Ringing the Randall-Sundrum braneworld: metastable gravity wave bound states 
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In the Randall-Sundrum scenario, our universe is a 4-dimensional 'brane' living in a 5-dimensional 
bulk spacetime. By studying the scattering ol bulk gravity waves, we show that this brane rings with 
a characteristic set of complex quasinormal frequencies, much like a black hole. To a bulk observer 
these modes are interpreted as metastable gravity wave bound states, while a brane observer views 
them as a discrete spectrum of decaying massive gravitons. Potential implications of these scattering 
resonances are discussed. 
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Introduction It is generally true that in many 
branches of physics, the best way to learn about an ob- 
ject is to hit it with something and analyze what comes 
back at you. Scattering experiments can either be real 
exercises done in the laboratory, such as in particle or 
condensed matter physics, or they can be of the virtual 
type, where the behaviour of radiation striking a target 
is simulated on a computer. The latter technique has 
proved to be an invaluable tool in the analysis of the 
fundamental structure of black holes, which are a class 
of targets not readily available in the lab. Numerical 
calculations of how gravity wave pulses bounce off black 
holes have revealed that they possess discrete quasinor- 
mal modes (QNMs) of vibration pj, which encode both 
important observational features for gravity wave detec- 
tors and valuable clues as to the canonical quantization 
of general relativity (GR) 0] ■ 

The success achieved for black holes prompts us to 
look for other virtual gravity wave scattering targets. 
A prime candidate comes from the Randall-Sundrum 
(RS) braneworld scenario This 5-dimensional phe- 
nomenological model is motivated by non-perturbative 
ideas from string theory, which postulate that we are 
living on a hypersurface (or 'brane') embedded within 
a higher-dimensional spacetime. One of the two vari- 
ants of the model involves one flat brane and an infinite 
fifth dimension. Usually, models with non-compact ex- 
tra dimensions are not viable because particles and fields 
exhibit higher-dimensional behaviour that is inconsistent 
with what we see in the world around us. However, the 
RS model escapes this fate via several key assumptions. 
First, standard model fields are confined on the brane, 
and hence behave in an entirely 4-dimensional manner. 
Second, although gravity is allowed to propagate in the 
bulk, at low energies it is effectively trapped on the brane. 
This is achieved by a negative cosmological constant in 
the bulk and enforcing the I2 reflection symmetry across 
the brane. The result is that the graviton has a con- 
tinuum of modes that have an effective 4-dimensional 
mass, and a single localized and normalizable massless 
mode. The latter 'zero' mode reproduces GR at low 
energies, while the former 'Kaluza-Klein' modes predict 
high-energy corrections to it. 

In this article, we consider the problem of scatter- 



ing gravity waves off the brane in the RS scenario. 1 We 
will demonstrate that just like black holes, the brane 
possesses a discrete set of QNMs that appear as scat- 
tering resonances; i.e., complex poles of the scatter- 
ing matrix. Since scattering resonances are essentially 
metastable bound states, this result has an important 
phenomenological implication for the RS one-brane sce- 
nario: namely, within the continuum of massive Kaluza- 
Klein graviton states, there are a discrete set of resonant 
modes that remain localized near the brane with definite 
lifetimes. After confirming the existence of these modes 
both analytically and numerically, we discuss how they 
might impact on various aspects of the braneworld sce- 
nario. 

Gravity wave master equation The RS scenario 
consists of an AdSs spacetime with a geometric defect 
(brane) at y — 0. The 5-metric is 

ds 2 = a 2 {y)^ v dx^dx v + dy 2 , a(y) = e^' 1 . (1) 

As is well known, in the RS gauge perturbations of the 
5-metric are written as 

9ab - 9ab + d a x^b\x v V, ^V-0 = c^V> (2) 

subject to the boundary condition dyh^ v — —2£~ 1 h^ l/ at 
y — 0. The simplest class of perturbations are the tensor 
modes, with harmonics — 1,2,3): 



V 2 T ( „ fc) - -fc 2 T„ (fe) 



d i T {k) = = TrT (fc) . 



I] I] 
Then, we express the metric perturbation as 



(3) 



Shij = 



d 3 k 
(2^r)3 



\y\l 2 ^ k {t,y)T^\ 5^ = 0. (4) 



Putting this ansatz into the linearized Einstein equations 
SG a b — (6/£ 2 )Sg a b yields that the $fc expansion coeffi- 
cients satisfy a master wave equation for y > 0: 

= [0 2 - 8 2 x + V k (x)] $ fc , V k (x) =k 2 + fx- 2 , (5a) 



= (d x 



$h at x = 1. 



(5b) 



1 We note that wave propagation (gravitational and otherwise) 
around 4-dimensional domain walls and supergravity 'D- 
brane' solutions 0] has been considered by several authors. 
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-38(x - 1) ^ 



the Schrodinger-like equation: uj 2 (f>k 



V k (x)(j) k . 



the boundary condition l|5b(l at the branc: 



If lu ^= k , this admits an exact solution that satisfies 



(7) 



H 



(1 V)i4 2 Vx) 



H 



(2 V)h^(f) 



Here, and H£ 2) arc Hankcl functions of the first and 
second kind respectively, and v is their order. The dimen- 
sionlcss /i parameter is defined via the dispersion relation 



FIG. 1: A schematic of the brane potential and incident grav- 
ity wave pulse in our scattering problem 

Here, we have defined dimensionless coordinates r = t/i 
and x — e v ^ > 1, as well as the dimensionless wavenum- 
ber k = k£. Note that the boundary condition l|5b|l and 
the x < 1 region can be incorporated into the potential 
via a delta-function and the imposition of mirror symme- 
try: 

V k {x) = k 2 + fx- 2 - 35(x- 1), x = \x-l\ + l. (6) 

So, the problem of perturbing the RS one-brane scenario 
is equivalent to the problem of wave propagation in a 
one-dimensional potential. 

Quasinormal resonances The situation we are in- 
terested in is shown in Figure ^ We regard the spatial 
wavenumber n of the perturbation as given, and we con- 
sider a compact pulse of right-moving gravitational radi- 
ation centered about xq > 1 at t = 0. This pulse will of 
course strike the brane and scatter, and we are interested 
in the late-time features of the reflected waveform. Now, 
since the potential governing gravity waves © is centered 
about an attractive delta- function, we can reasonably ex- 
pect some of the incident wave to become trapped and 
remain localized near the brane forever, while the rest of 
energy will bounce off and propagate to infinity. 

However, all the reflected energy need not be promptly 
reflected away from the brane. A well-known phe- 
nomenon in spectral theory is that of a scattering res- 
onance [||. This occurs when certain frequency com- 
ponents of the initial pulse become temporarily trapped 
within the potential. One can think of these modes as 
metastable bound states: they are initially localized, but 
they eventually leak out to infinity. In black hole per- 
turbation theory, these metastable states are precisely 
the QNMs. In that problem, the master equation has 
the form of l|5a|) with a potential consisting of a smooth 
barrier centered about the photon-sphere. The late time 
behaviour of a pulse after bouncing off this barrier is 
dominated by QNMs with exp(ito n T) time-dependence, 
where {w n } is a discrete set of complex frequencies with 
Im uj n > 0; i.e., quasinormal oscillations are damped in 
time. 

The question is now: does the gravitational wave po- 
tential in the RS wave equation (|5a|l support QNMs? To 
answer, we first assume $fc = e lu}T 4>k(x), which yields 



^ 2 =u 2 -k 2 , |Arg(//)| < 7T. (8) 

Conventionally, \x is assumed to be real and positive, and 
is related to the conventional Kaluza-Klein mass of the 
fluctuation via fi = ml. However, the key to finding a 
resonant mode is to let /j, range over the complex plane. 
Consider the asymptotic expansion of the Hankel func- 
tions: 

<f>k(x) cx B^\a)e-^ x - Hf ) (ji)e +i " a> as x — oo. (9) 

Depending on the sign of Re /i relative to the sign of 
Re u>, each term corresponds to either a left or right 
moving wave at infinity. The next step in the QNM pro- 
gramme is to demand that the asymptotic form of 4>k is 
'purely outgoing' [J, which in this context means that 
we need to set the coefficient of the left-moving wave to 
zero. 2 Formally, we are searching for a complex pole of 
the scattering matrix Sk(w), which is defined as the ratio 
of the amplitudes of outgoing and ingoing radiation at 
infinity. Clearly, for a pole we need fi to be a zero of 
either Hj 1 ^ or H^ 2 -* . For reference, we quote the first four 
zeros of H^: 

Ht = -0.419 -0.577 i, /x 2 = -3.832 - 0.355 i, 
fi 3 = -7.016 - 0.350 i, /z 4 = -10.174 - 0.348 i. 

(2) 

The corresponding zeros of H^ are simply the complex 
conjugates p, n , and it is easy to confirm that all zeros of 
both Hankel functions have negative real parts. There- 
fore, if Re oj > we must set H^ [a) = in order to have 
a purely outgoing wave; that is, we choose fi £ {l^n}- 

(2) 

Conversely, if Re ui < we need H^ [u) = 0, or in other 
words [i S {/]„}. From this, it follows that there are two 
families of complementary quasinormal frequencies: 

^(k) = +y/l% + HfV„) = 0, (10a) 

4 2 >( K ) = Hf^Ho. ( iob) 



2 The fact that QNMs are defined by a boundary condition 'at in- 
finity' tells us that they are not relevant to two-brane scenarios 
where the extra dimension is effectively finite. Generally speak- 
ing, in order for a system to exhibit metastable resonances there 
must be some type of mechanism to remove energy from the 
system, which only really happens in the one-brane case. 
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FIG. 2: The first twenty quasinormal frequencies of the brane 
as the spatial wavenumber n is varied from (red) to 25 (pur- 
ple). Those in the right half of the complex plane correspond 
to lj! 1 ' , while those in the left half are w„ 



The variation of the first twenty QNMs with k is depicted 
in Figure El Note that in order to evaluate the square 
root, we need to take the same branch cut as that of the 
Hankel function; i.e., the negative real axis. In particular, 
this means \fi^ = — u n . 

Before moving on, we should mention that the Randall- 
Sundrum model is not the only 5-dimensional scenario 
that exhibits resonant behaviour. The Dvali-Gabadadze- 
Porrati (DGP) model features a Minkowski bulk sur- 
rounding a brane with 'induced gravity' pj. There is 
no zero-mode in this model, but the massive modes sum 
together in such a way to approximately recover New- 
ton's law on small scales. This superposition of mas- 
sive modes is interpreted as a metastable 4-dimensional 
graviton that decays into the bulk after a finite streaming 
length along the brane. To the best of our knowledge, the 
QNM spectrum of the DGP model has not been specifi- 
cally investigated, although for certain branch cut choice 
one can find DGP resonances with Imw < 0; i.e., unsta- 
ble fluctuations |g. An interesting and perhaps fruitful 
exercise would be to repeat the current calculation for 
that model. 

Numerical results So far, we have concentrated on 
modes that are outgoing at spatial infinity. But the RS 
potential also supports a normalizable bound state: 



(11) 



The existence of this so-called 'zero-mode' is what gives 
the RS model its main appeal, since it gives rise to GR 
at low energies. How does the zero mode manifest itself 
in our scattering problem? To answer, we first look at a 
numeric simulation of a Gaussian pulse incident on the 
brane in Figure |3J This is obtained by numerically inte- 
grating (|5a|l for x > 1 and with the boundary condition 
(|5b|) . Our assumed initial data is at x — 10 and has a 
full-width at half maximum (FWHM) of 5, and we take 
k = 7r/10. We see that after the pulse hits the brane, 
part of the energy is reflected and propagates off to in- 
finity while the other part remains trapped on the brane 
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FIG. 3: Numeric waveform for « = 7r/10 and an incident 
Gaussian pulse centered at x = 10 with a FWHM of 5 



at x = 1. Indeed, a close examination of the brane per- 
turbation shows that after the initial pulse, the waveform 
is an undamped sinusoid with frequency 7r/10. On the 
other hand, the waveform at x = 100 shows that the re- 
flected pulse damps to zero as r — > oo. In other words, 
the incident pulse excites a perpetual and spatially lo- 
calized perturbation on the brane; i.e., the zero- mode. 

The one thing missing from the situation depicted in 
Figure |21 is any hint of resonant phenomena. The Gaus- 
sian pulse immediately excites the zero mode when it 
strikes the brane, which suggests any quasinormal reso- 
nances are sub-dominant to the undamped bound state. 
Some experimentation reveals that this is true for many 
different types of Gaussians and values of k. The reason 
is simple: in Fourier space a Gaussian pulse always has a 
significant low frequency component that will efficiently 
excite the zero mode. In order to see the QNM influ- 
ence directly, we should really consider coherent plane 
wave initial data with a frequency high enough to pref- 
erentially excite the first QNM. But we cannot model 
an infinite plane wave numerically, so we instead use a 
finite wavetrain, which is essentially a segment of the left- 



moving wave sin(ix)o' r + IM)x), with u>q 



Mo 



If such 



a wavetrain is N cycles long, its spectral width is on the 
order of Aoj = ujq/N. It is also useful to look at the 
reflected wave from an off-brane vantage point, because 
the influence of the zero mode is strongest at x = 1. 
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In Figure we show the perturbation evaluated at 
x = 10 for this type of initial data with k = tt/10 and 
u>o = 0.6. The lefthand edge of the wavetrain is initially 
at x = 10, and it is 10 wavelengths long. For this choice of 
k, the first quasinormal mode is at uj^ = 0.463 + 0.523 i. 
In this case, the gravitational waveform is analogous to 
the response of a damped mechanical oscillator subject 
to a driving force. The natural frequency of the sys- 
tem is that of the zero mode u> — 7r/10 ~ 0.314, so the 
wavetrain represents a high frequency forcing term. As 
the wavetrain reflects off the brane, the perturbation is 
obliged to oscillate at uj = 0.6. By t ~ 150, the entire 
reflected wavetrain has travelled past our vantage point, 
and the driving force is effectively 'switched off.' But the 
waveform does not immediately revert to zero mode vi- 
brations, there is instead a 'ringdown phase' where the 
brane radiates away the high frequency components of 
the driving force. The upper inset of Figure^ shows the 
amplitude squared of the waveforms's Fast Fourier Trans- 
form (FFT) during this ringdown. The power spectrum 
is peaked at the real part of the first QNM at 0.463, 
which is fairly clear evidence of quasinormal ringing of 
the brane. (Note that the damping of a QNM means 
that its power spectrum has a significant width.) Af- 
ter the ringdown phase, the FFT (not shown) is sharply 
peaked about 7r/10, suggesting strong domination of the 
zero mode. In the standard lexicon of black hole per- 
turbation theory, the part of the waveform that follows 
quasinormal ringing is known as the 'wave-tail.' So in the 
RS braneworld, the wave tail essentially consists of zero 
mode oscillations. The lower inset in Figure0]shows that 
the amplitude of the zero mode signal becomes constant 
only for late times. Note that a slow variation in the 
amplitude of oscillatory wave-tails is a familiar feature of 
potentials that are non- vanishing at infinity 0- In the 
future, we will report on the Randall-Sundrum wave-tail 
in detail. 

Implications As mentioned above, the principle im- 
plication of the existence of QNMs in the RS onc- 
brane scenario is the prediction of a discrete set of 
metastable massive modes within the continuous Kaluza- 
Klein tower. 3 If I = 0.1 mm, which is an upper value set 
by experimental limits on deviations from Newton's law 
|12| . the lowest order resonance has mass 8.386 x 10 -4 eV 
and a half life of 4.002 x 10~ n seconds; i.e., it is light and 
short-lived. Ideally, one would look for such resonances 
in particle accelerators, but in the RS one-brane case the 
coupling between matter and the Kaluza-Klein states is 
quite small since the 5-dimensional Planck mass is > 10 s 
TeV. This scale suggests that the best place to look for 
QNMs is in the high-energy regime associated with the 
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3 Alternatively, one could view the resonant modes as spin-2 fields 
obeying a complex- valued dispersion relation; i.e., travelling in 
an absorptive medium. Details of such an interpretation can be 
found in Ref. 
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FIG. 4: Numeric waveform evaluated at x = 10 for an incident 
wavetrain with u>o — 0.6 and 10 wavelengths long. We take 
K = tt/10. The upper inset shows the power spectrum of the 
waveform during the ringdown phase, while the lower inset 
shows its late time behaviour 



early universe; in particular, one can expect the stochas- 
tic gravity wave background [10j to carry potentially ob- 
servable echoes of metastable masses. Of course, the ex- 
pansion of the brane universe at that epoch will certainly 
change the spectrum of quasinormal resonances. How- 
ever, it is reasonable to assume that very high-frequency 
QNMs will 'see' a moving brane to be stationary, and 
hence our results from the static case will be applicable. 

It is interesting, and perhaps instructive, to see what 
contribution the quasinormal modes make to the effec- 
tive braneworld Einstein equations. Using the Shiromizu- 
Maeda-Sasaki formalism [lg, we see that QNMs induce 
'dark radiation' on the brane essentially the pullback 
of the 5-dimensional Weyl tensor onto the 4-dimensional 
spacetime. The QNM contribution to the stress energy 
of this effective matter field is 



rrrp(fe) 
R /J.V ' 



(12) 



where A n are constants describing the degree to which 
mode is excited. Note the explicit dependance on the 
complex mass /i„; the massless zero-mode does not enter 
into the dark radiation. Hence, the quasinormal reso- 
nances can be interpreted as a discrete set of brane mat- 
ter fields prone to decay by tunnelling into the bulk. Re- 
turning to the cosmological case, this implies that the 
principal high frequency bulk effects on early universe 
cosmological perturbations are via 'matter' contributions 
of the form 112(1 . The quantitative effects of these fields 
on the cosmic microwave background and large scale 
structure are yet to be calculated. 

Finally, Eq. (|12|l hints at a possible holographic in- 
terpretation of these resonances. According to the 
AdS/CFT conjecture, the Randall-Sundrum model is 
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dual to a large- N conformal field theory (CFT) with cut- 
off living on the brane 01 • Specifically, the retarded bulk 
graviton Green's function is identified with the generat- 
ing functional of the CFT. Now, the metastable reso- 
nances found here are complex poles of the Green's func- 
tion, and hence will play an important role in the bound- 
ary theory. In particular, note that £^ v in the gravita- 
tional theory is precisely the stress-energy tensor of the 
CFT, which suggests that the gravitational QNMs de- 
scribe how the boundary theory relaxes to equilibrium 
after some initial disturbance, and also give the timescale 
for such a relaxation. 

Summary Just like a black hole, the RS braneworld 
'rings' with a discrete set of complex frequencies. These 
metastable bound states, which look like decaying mas- 



sive gravitons to 4-dimensional observers, should play an 
important role in any brane phenomena involving bulk 
degrees of freedom. An interesting question is: how do 
QNMs manifest themselves in generalizations of the RS 
scenario? We have recently reported on the resonant be- 
haviour of a spherical Einstein static brane encircling a 
Schwarzschild- AdSs black hole > but the key challenge 
is characterizing the quasinormal excitations of a moving 
cosmological brane. We hope to investigate this in the 
future. 
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